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ABSTRACT 

We test models for the evolution of neutron star (NS) magnetic fields (B). 
Our model for the evolution of the NS spin is taken from an analysis of pulsar 
timing noise presented by Hobbs et al. (2010). We first test the standard model 
of a pulsar's magnetosphere in which B does not change with time and magnetic 
dipole radiation is assumed to dominate the pulsar's spin-down. We find this 
model fails to predict both the magnitudes and signs of the second derivatives 
of the spin frequencies (u). We then construct a phenomenological model of the 
evolution of B, which contains a long term decay (LTD) modulated by short term 
oscillations (STO); a pulsar's spin is thus modified by its B-evolution. We find 
that an exponential LTD is not favored by the observed statistical properties of 
v for young pulsars and fails to explain the fact that v is negative for roughly 
half of the old pulsars. A simple power-law LTD can explain all the observed 
statistical properties of i). Finally we discuss some physical implications of our 
results to models of the S-decay of NSs and suggest reliable determination of 
the true ages of many young NSs is needed, in order to constrain further the 
physical mechanisms of their I?-decay. Our model can be further tested with 
the measured evolutions of z> and i) for an individual pulsar; the decay index, 
oscillation amplitude and period can also be determined this way for the pulsar. 

Subject headings: stars: neutron — pulsars: general — magnetic fields — methods: 
statistical 



1. Introduction 

Many studies on the possible magnetic field decay of neutron stars (NSs) have been 
done previously, based on the observed statistics of their periods (P) and period derivatives 



1 National Astronomical Observatories, Chinese Academy of Sciences, Beijing, 100012, China; 
zhangsn@ihep .ac.cn 

2 Key Laboratory of Particle Astrophysics, Institute of High Energy Physics, Chinese Academy of Sciences, 
Beijing 100049, China 



- 2 - 



(P). Some of these studies relied on a population synthesis of pulsars (Holt & Ramaty 1970; 
Bhattacharya et al. 1992; Han 1997; Regimbau & de Freitas Pacheco 2001; Tauris & Konar 
2001; Gonthier et al. 2002; Guseinov et al. 2004; Aguilera et al. 2008; Popov et al. 2010); 
however, no firm conclusion can be drawn on if and how the magnetic fields of NSs decay 
(e.g., Harding & Lai 2006; Ridley & Lorimer 2010; Lorimer 2011). Alternatively some other 
studies used the spin-down or characteristics ages, r c = (P — P )/2P (P is the initial spin 
period of the given pulsar), or r c = P/2P if P 3> Pq, as indicators of the true ages of NSs and 
found evidence for their dipole magnetic field decay (Pacini 1969; Ostriker & Gunn 1969; 
Gunn & Ostriker 1970). However, as we have shown recently (Zhang & Xie 2011), their 
spin-down ages are normally significantly larger than the ages of the supernova remnants 
physically associated with them, which in principle should be the unbiased age indicators of 
these NSs (Lyne 1975; Geppert et al., 1999; Ruderman 2005). It has been shown that this 
age mismatch can be understood if the magnetic fields of these NSs decay significantly over 
their life times, since the magnetic field decay alters the spin-down rate of a NS significantly, 
and thus cause the observed age mismatch (Lyne 1975; Geppert et al. 1999; Ruderman 2005; 
Zhang & Xie 2011). 

Pulsars are generally very stable natural clocks with observed steady pulses. However 
significant timing irregularities, i.e., unpredicted times of arrival of pulses, exist for most 
pulsars studied so far (see Hobbs et al. 2010, for an extensive reviews of many previous 
studies on timing irregularities of pulsars). The timing irregularities of the first type are 
"glitches", i.e., sudden increases in spin rate followed by a period of relaxation; it has been 
found that the timing irregularities of young pulsars with r c < 10 5 years are dominated 
by the recovery from previous glitch events (Hobbs et al. 2010). In many cases the NS 
recovers to the spin rate prior to the glitch event and thus the glitch event can be removed 
from the data satisfactorily without causing significant residuals over model predictions. 
However in some cases, glitches can cause permanent changes to both P and P of the NS, 
which cannot always be removed from the data satisfactorily, for instance, if the glitch has 
occurred before the start of the data taking and/or long term recoveries from glitch events 
cannot be modeled satisfactorily. These changes can be modeled as a permanent increase of 
the surface dipole magnetic field of the NS; consequently some of these NSs may grow their 
surface magnetic field strength gradually this way and eventually have surface magnetic field 
strength comparable to that of magnetars over a life time of 10 5-6 years (Lin & Zhang 2004). 
These are the first studies that linked the timing irregularities of pulsars with the long term 
evolution of magnetic fields of NSs. 

The timing irregularities of the second type have time scales of years to decades and 
thus are normally referred to as timing noise (Hobbs et al. 2004; 2010). Hobbs et al. (2004, 
2010) carried out so far the most extensive study of the long term timing irregularities of 366 
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pulsars. Besides ruling out some timing noise models in terms of observational imperfections, 
random walks, and planetary companions, Hobbs et al. (2004, 2010) concluded that the 
observed v values for the majority of pulsars are not caused by magnetic dipole radiation or 
by any other systematic loss of rotational energy, but are dominated by the amount of timing 
noise present in the residuals and the data span. Some of their main results that are modeled 
and interpreted in this work are: (1) All young pulsars have v > 0; (2) Approximately half 
of the older pulsars have v > and the other half have v < 0; (3) Quasi-oscillatory timing 
residuals are observed for many pulsars; and (4) The value of v measured depends upon the 
data span and the exact choice of data processed. 

Proper understanding of these important results can lead to better understanding of 
the internal structures and the magnetospheres of NSs. Since the spin evolution of pulsars 
cannot be caused by any systematic loss of rotational energy (Hobbs et al. 2004; 2010), it is 
natural to introduce some oscillatory parameters. Indeed, some of the observed periodicities 
have been suggested as being caused by free-precession of the NS (Stairs et al. 2000). 
Some strong quasi-periodicities have been identified recently as due to abrupt changes of 
the magnetospheric regulation of NSs, perhaps due to varied particle emissions (Lyne et 
al. 2010). It was suggested that such varied magnetospheric particle emissions is also 
responsible for the observed long term timing noise (Liu et al. 2011). Alternatively it has 
been suggested that the propagation of Tkachenko waves in a NS may modulate its moment 
of inertia, and thus produce the observed periodic or quasi-periodic timing noise (Tkachenko 
1966; Ruderman 1970; Haskell 2011). Because all of the above processes can only cause 
periodic or quasi-period modulations to the observed spin of NSs, all these modulations can 
be modeled phenomenologically as some sort of oscillations of the observed surface magnetic 
field strengths of NSs. It is thus reasonable to attribute the observed wide-spread quasi- 
oscillatory timing residuals to the oscillations of the magnetic fields of NSs. As shown in this 
work, the oscillating magnetic fields can also naturally explain the fact that approximately 
half of the older pulsars have v > and the other half have v < 0, when the long term 
magnetic field decay does not dominate the spin evolutions of pulsars. This is the main 
advantage of this model, and thus justifies invoking oscillating magnetic fields for modeling 
the long term spin evolutions of pulsars. 

In this work we attempt to use the known statistical properties of the spin evolution of 
the pulsars to test some phenomenological models for the evolution of the dipole magnetic 
field of a NS. In sections 2-5, we test the following three models: (1) the standard magnetic 
dipole radiation model with constant magnetic field, (2)** long-term exponential decay with 
an oscillatory component and (3) long-term power-law decay modulated with oscillations. 
We find the last model is the only one compatible with the data, with the power-law index 
a > 0.5 and the oscillation amplitude of around 5 x 10~ 5 for an oscillation period of several 
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years. Finally we discuss the physical implication of our results, summarize our results and 
make conclusions and discussions, including predictions of the model for further studies. All 
timing data of pulsars used in this work are taken from Hobbs et al. (2010). 



2. Testing the standard magnetic dipole radiation model 

Assuming pure magnetic dipole radiation as the braking mechanism for a pulsar's spin- 
down, we have, 

inn = -^rV (i) 

or, 

v = -AB 2 v 3 , (2) 

where A = - 6e3/ J is normally assumed to be a constant, B is its dipole magnetic field at its 
magnetic pole, R is its radius, I is its moment of inertia. Assuming B = 0, we have, 

v = 3v 2 /v. (3) 

This model thus predicts v > 0, against the fact that many pulsars show i) < 0. This is 
a major failure of this model, or any other models with monotonic spin-down, as pointed 
out in Hobbs et al. (2010). In Fig. CEJ we plot the comparison between the observed i) and 
that predicted by Eq. ([3]). It is clear that this model under-predicts the amplitudes of i) by 
several orders of magnitudes for most pulsars. 

Assuming B ^ 0, we have, 

v = 3v 2 /v + 2vB/B. (4) 

From the above analysis, it is obvious that the second term in Eq. (T4|) must dominate the 
observed i) (i.e., 3v 2 jv <C \2vB/B\), for both its magnitude and sign. Define K = AB 2 , we 
have, 

We define the evolution time scale of K, tk = K/\K\ = B/2\B\, we have, 

i ™ I (R\ 

tk = \— 2 — • (6) 

Fig. [2] shows the observed tk as a function of r c . Clearly tk -C t c for most pulsars. This 
suggests that if the measured i) values are caused by magnetic field evolution, then the 
evolution must be significant for most pulsars in the sample. Fig. |2] also shows that for all 
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Fig. 1. — Observed correlation between v and 3u 2 /u. The prediction of the standard mag- 
netic dipole radiation model, i.e., v = 3z> 2 /V is shown as the dashed line, which under-predicts 
significantly the magnitudes of v for most pulsars and also cannot explain i) < for nearly 
half of the pulsars. 
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Fig. 2. — Correlation between the time scale of K variation, tk = K/\K\, and the character- 
istic age, r c , of pulsars; r K = r c is shown as the dashed line. All pulsars show clear variation 
of K. 



-6 - 



young pulsars with r c < 10 5 yr, K < 0, i.e., their magnetic moment field decays with time. 
For the rest of pulsars, there are almost equal numbers of pulsars with positive or negative 
K. If one assumes that most pulsars are born with some common properties and they then 
evolve to have the observed statistical properties, it is then necessary that at least some 
pulsars should have both K < and K > at different times, i.e., K oscillates. 

We therefore model phenomenologically the evolution of the observed dipole magnetic 
field of a pulsar by a decay modulated by the dominated component of oscillations, 

B(t) = B d (t)(l + ksm(<f> + 27T^)), (7) 

where t is the pulsar's age, k, (ft and T are the amplitude, phase and period of the oscillating 
magnetic field, respectively. Sometimes it might be necessary to take multiple oscillation 
components since multiple peaks are often seen in the power spectra of the timing residuals 
of many pulsars (Hobbs et al. 2010). Assuming that all the oscillation time scale is much 
less than the decay time scale (as will be discussed and justified later), i.e., T <C B d /B d , we 
have, 

B~B d + B d fcos(0 + 2n^), (8) 

where / = 2nk/T. For simplicity here we only consider the amplitude of the oscillation term 
to model the statistical properties of the timing noise of pulsars, i.e., 

B = B d ± fB d . (9) 

Since B < can produce v > 0, a viable model should at least have \B d \/B d > f for young 
pulsars and vice versa for old pulsars, in order to produce positive v for young pulsars, but 
either positive or negative v for old pulsars with almost equal chances, as discussed above. 



3. Testing Exponential Decay Model with Young Pulsars 

Assuming B d = B exp(— t/r d ), Eq. (jUJ) becomes, 

B = B d (-l/r d ±f). (10) 

Substituting Eq. ffTUl) into Eq. (TJ]) and ignoring its first term (since 3z> 2 /z/ \2i/B/B\), we 
have, 

i>~-2u(l/r d ±f). (11) 

Since v > for most young pulsars, we assume l/r d > f, i.e., v ~ —2vjT d , or, r d ~ —2vjv 
for young pulsars. In Fig. [31 we show the observed v and the inferred r d as a function of 
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the observed u, for young pulsars with r c < 10 6 yr and i) > 0. Clearly a single value or 
distribution of r d cannot describe the data adequately. Another difficulty of this model is 
that i) < for the nearly half of the pulsars cannot be explained with the same assumption 
of l/r d > /. We thus conclude that a simple exponential decay model is not favored by 
data. 



4. Testing Power-law Decay Model with Young Pulsars 

Assuming B d = B (t/to)~ a , Eq. (jUJ) becomes, 

B = B d (-a/t±f). (12) 

Substituting Eq. ffl2l into Eq. (j3j) and ignoring its first term (since 3z> 2 /z/ <C \2i>B/B\), we 
have, 

D~-2v(a/t±f), (13) 

where t is the real age of the pulsar. Without other information about t, here we replace it 
with the magnetic age of a pulsar t = t (B / B)^ 1 ^^; we then have, 

i>~ v (-v) l+l3 /is 3 P ±2vf, (14) 

where = l/2a, r) = (3.3 x 10 19 / B ) 2 ^ / fit and B = 3.3 x 10 19 y/p? G is assumed. 
Since v > for most young pulsars, we assume a/t > f, i.e., 

v ~ r](-z>) 1+/3 /z/ 3/3 , (15) 
for young pulsars. For the i-th pulsar, we have, 

Vi = Vi/i-isY+Pu-W, (16) 
We then calculate fit the root mean squares (rms) for each value of ft = 1 /2a, 

rms = — Jy^iVi-Vm) 2 , (17) 

where n is the number of pulsars and r] m is the median value of all 77,. In Fig. HJ^a) we 
show rms as a function of a, indicating that a > 0.5 is required, but its upper limit is 
unconstrained. In Fig. HJ^a, b) we also show the ranges of r/ m and the corresponding to and 
B . 

The uncertainty in determining the exact value of a is due to the unknown real ages 
of these pulsars; using the magnetic ages as their ages introduces a strong coupling between 



-8 - 



the parameters in the magnetic field decay model and their ages. This is manifested in the 
correlations of a with rms, r/ m , t Q and B Q . Nevertheless, their ranges are reasonable with 
a > 0.5. In Fig. [S^a, b), we show the comparison between the predictions of Eq. f[T5"j) and 
data with a = 0.5 and a = 1.0, respectively. In both cases the model can describe the data 
reasonably well. We thus conclude that a simple power-law decay model is favored by data. 

5. Testing the Model of Power-law Decay Modulated by Oscillations with All 

Pulsars 

Eq. (12 ip shows that the oscillation term / dominates both the signs and magnitudes of 
v for old pulsars with t > a/ f; this explains naturally why there are almost equal numbers 
of pulsars with v > or v < 0. This agrees with the observation that a single old pulsar 
can be observed to show either v > or v < 0, depending on which segment of data is used 
(Hobbs et al. 2010). We also confirmed this with Monte-Carlo simulations in a preliminary 
study (Figs. 5 and 6 in Zhang & Xie 2012); a more detailed study will be presented elsewhere 
(Xie & Zhang 2012). Fig. S(a) shows that rms remains almost constant with a > 1; we thus 
restrict our analysis in the following for only two cases: a = 0.5 and a = 1.0. 

Eq. ( 12~TI) can rewritten as, 

i) v W ~ v (-v) 1+f> ± 2z>z/ 3/3 /, (18) 

suggesting that / ^ causes deviations to the expected correlation between i)u 3f3 and z>, as 
shown in Fig. [HI Choosing the same range of /, the data for both positive and negative v can 
be explained; r\ is taken as rj m determined from young pulsars with i) > and r c < 2 x 10 6 yr, 
for a = 0.5 and a = 1.0, respectively. For each pulsar in the sample, / can also be determined 
from Eq. (l2~Tj) or Eq. ( TT8|) . as shown in Fig. [7J clearly / is not correlated with either v or v 
(we also find that / is not correlated with either v or r c , which is not shown here for brevity) . 
This indicates that the statistical properties of the timing noise of all kinds of pulsars can 
be explained with the same oscillation model of similar oscillation amplitudes. 

In Fig. [HI the distributions of f /a, for all pulsars with positive or negative i), are shown 
in the cases of a = 0.5 and a = 1.0, respectively. All distributions show well-defined peaks 
and consistent to each other and the peak value of f/a is around 3 x 10~ 12 ; we also find 
that choosing different values of a with 0.5 < a < 2.0 always gives similar f/a. Finally in 
Fig. [8] we show the comparisons between the observed v and that predicted by Eq. ([2~T]) with 
f/a = 3 x 10~ 12 for a = 0.5 and a = 1.0, respectively. This simple model can describe the 
general properties of the data over about nine orders of magnitudes of v. 
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6. Physical Implications 
6.1. The braking index 

Eq. dSJ) can be rewritten as, 

K=- 3 (--l), (19) 

where rib = vvjv 2 is the conventionally defined braking index. Eq. (119]) means that K > 
when nb < 3, including < when i> < 0. This agrees with many previous works reporting 
that nb < 3 indicates magnetic field increase, and vice versa (e.g., Blandford & Romani 
1988; Chanmugam & Sang 1989; Manchester & Peterson 1989; Lyne et al. 1993; Lyne et al. 
1996; Johnston & Galloway 1999; Wu et al. 2003; Lyne 2004; Chen & Li 2006; Livingstone 
et al. 2006; Middleditch et al. 2006; Livingstone et al. 2007; Yue et al. 2007; Espinoza et 
al. 2011). However, these works have reported very few cases that rib is measured reliably. 

In Fig. [TUl we show the distribution of in this sample. There is no clustering around 
nb = 3 at all, and the values positive rib and the absolute values of negative rib exhibit almost 
identical distributions. The numbers of positive or negative values of rib are 193 and 173, 
respectively. Similarly the numbers of rib > 3 or rib < 3 (including all negative rib) are 192 
and 174, respectively. Statistically these numbers are indistinguishable. Therefore about 
half of pulsars in this sample are observed with increasing or decreasing magnetic fields. 



6.2. Magnetic field decay 

Several physical mechanisms have been proposed for magnetic field decay in NSs, e.g., 
ohmic dissipation, Hall effect, and ambipolar diffusion. We can therefore have (Goldreich & 
Reisenegger 1992; Heyl & Kulkarni 1998), 

dB _ 1 + 1 + 1 

dt ^ohmic ^ambip ^hall 

where t hmic ~ 2 x 10 11 ^|(^ £ -) 3 yr, ^mbip ~ 3 x 10 9 -^^- yr (for the solenoidal component), 
^ambip ~ ^5^2 (l + 5 x lO^L^T!) yr (for the irrotational component), and than ~ 5 x 

10 ~bI7"( _£!_ ) y r ' wnere B\2 is the surface magnetic field strength in units of 10 G, Tg 
denotes the core temperature in units of 10 8 K, p nuc = 2.8 x 10 14 g cm -3 is the nuclear 
density, and L 5 is a characteristic length scale of the flux loops through the outer core in 
units of 10 5 cm. 

Assuming T 8 does not vary or only varies slowly with time (Pons et al. 2009), the ohmic 
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dissipation will produce an exponential decay of the magnetic field, B = BQe~ l / Td . The other 
two mechanisms all produce power-law decay, B = B (t/t )~ a , where a = 1 or a = 0.5 
for the Hall effect or ambipolar diffusion (for t <C 10 12 yr), respectively. We recently found 
that this model with constant core temperatures of NSs can explain well the observed age 
mismatch of NSs discussed above, and the ambipolar diffusion effect seems to be dominant, 
though the statistics is rather limited (Zhang & Xie 2011; Xie & Zhang 2011) 

As shown in Fig. 03) and discussed in Sect. 3, the simple exponential decay cannot 
describe the data for young pulsars with r c < 10 6 yr. Therefore a constant core temperature 
(T§ = const) ohmic dissipation model for all magnetic field decay is ruled out. However, 
for a cooling NS, i.e., Tg decreases with time, the ohmic dissipation model would result in 
an increasing t$ with time, in rough agreement with that shown in Fig. ©(b). Indeed, 
an increasing t$ with time would also make the first term in Eq. ( TTTj) less important or 
negligible for old pulsars, thus explaining why there are almost equal number of pulsars with 
either positive or negative i). Of course an increasing t& with time may be able to mimic a 
power-law decay mathematically. 

As shown in Figs. (JH1H]) and discussed in Sects. 4 and 5, a power-law decay with 
a > 0.5 can explain all the observed statistical properties of the timing noise of pulsars. 
Unfortunately the exact value or even the upper limit of a cannot be constrained by the 
data, though a ~ 1 is preferred. For both the Hall effect and the ambipolar diffusion (for 
t <C 10 12 yr), a cooling NS implies faster magnetic field decay, and thus increasing a with 
time. 

We therefore conclude that the current data do not conflict the above discussed model 
of magnetic field decay of NSs, but also cannot constrain precisely how these effects regulate 
their magnetic field decay. 

6.3. Magnetic field oscillations 

As demonstrated in Sect. 5, invoking an oscillation component with //a ~ 3x 10~ 12 s -1 , 
both the amplitude and signs of i) can be well reproduced for all pulsars. Quasi-periodic 
oscillations with time scales of the order of years in the timing residuals of pulsars are found 
to be quite common (Hobbs et al. 2010). As we have shown recently, such quasi-periodic 
oscillations can be re-produced by an oscillating dipole magnetic field of a NS with an 
oscillation period of several years (Zhang & Xie 2012). Taking T = 3 years and a = 1, we 
have the oscillation amplitude of the observed magnetic field fc~5x 10~ 5 . 

Assuming that these oscillations are caused by precessing of the NS and the angle 
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between its magnetic and spin axes is 60 degrees, k ~ 5 x 10~ 5 requires a precession angle 
of about 0.003 degree, which is much smaller than the inferred 0.3 degree precession angle 
for PSR B1828— 11 from its periodic timing properties (Stairs et al. 2000). Therefore 
PSR B1828— 11 is probably an extreme example showing observed oscillating magnetic fields, 
though our result suggests that the same phenomenon exists in essentially all pulsars. Of 
course the apparent oscillations of their observed magnetic fields may also be caused by 
periodic or quasi-periodic changes of their magnetospheres (Lyne et al. 2010; Liu et al. 
2010), or even their moments of inertia caused by, for example, the Tkachenko waves in NSs. 
Regardless which mechanism is responsible for the oscillations, the mechanism appear to work 
in the same way for all types of NSs, since the distribution of f /a is quite narrow (especially 
when a = 1), and / does not appear to be correlated with any observed parameters of 
pulsars. 



7. Summary, conclusion and discussion 

In this work we have tested models of magnetic field evolution of NSs with the statistical 
properties of their timing noise reported by Hobbs et al. (2010). In all models the magnetic 
dipole radiation is assumed to dominate the spin-down of pulsars; therefore different models 
of their magnetic field evolution lead to different properties of their spin-down. Our main 
results in this work are summarized as follows: 

(1) We tested the standard model of pulsar's magnetosphere with constant magnetic 
field. We found that this model under-predicts the second derivatives of their spin frequencies 
(i>) by several orders of magnitudes for most pulsars (Fig. ([T])). This model also fails to 
predict i) < 0, which is observed in nearly half of the old pulsars. The standard model of 
pulsar's magnetosphere with constant magnetic field is thus ruled out, and v is found to be 
dominated by the varying magnetic field. 

(2) Statistically about half of the pulsars in this sample show increasing or decreasing 
magnetic fields. The magnetic field strength of most pulsars are found to either increase or 
decrease with time scales significantly shorter than their characteristic ages (Fig. fl2])). 

(3) We constructed a phenomenological model of the evolution of the magnetic fields of 
NSs, which is made of a long term decay modulated by short term oscillations. 

(4) A simple exponential decay is not favored by the apparent correlation of the inferred 
decay constant with their characteristic ages for young pulsars (Fig. ©(b)) and fails to 
explain the fact that i) is negative for roughly half of old pulsars. 
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(5) A simple power-law decay, however, can explain all the observed statistical properties 
of v (Figs. ()5] and [6])) and infer a narrow distribution of an oscillation amplitude of about 
5 x with a preferred decay index of about unity (Figs. ([7] and [5])). 

We thus conclude that long term magnetic field and short term oscillations are ubiq- 
uitous in all NSs and that the data support the model consisting of a power-law decay 
modulated by short term oscillations for the evolution of their magnetic fields. However we 
can only constrain the decay index to be a > 0.5 (Figs. fll]) and ©(b)). 

Our results are consistent with the models of magnetic field decay we discussed, in 
which ohmic dissipation, Hall effect, and ambipolar diffusion are possible effects causing 
magnetic field decay. However, the uncertainty in the decay index a makes it difficult to 
constrain further or distinguish between these effects. The problem in determining a is the 
unavailability of the true ages of pulsars. With the true age t available for some pulsars, a 
can be determined directly with Eq. (fl3|) ; / = can be safely assumed for young pulsars. 
Unfortunately only several pulsars in the sample of Hobbs et al (2010) have the ages of their 
supernova remnants available, thus preventing a statistically meaningful determination of a 
this way. We thus suggest that reliable determination of the true ages of a large sample of 
young NSs is needed, in order to constrain further the physical mechanisms of the magnetic 
field decay of NSs. This can also provide an independent test of our model. 

The main prediction of this model is that the spin evolution of each pulsar should follow, 

i>(t) ~ 2u(t)(~ + J2f* cos(& + 2tt|-)), (21) 

where both i)(t) and i>(t) should be time-resolved, rather than the averaged value conven- 
tionally reported over the whole data span. Hobbs et al. (2010) have shown that both the 
value and sign of the averaged v depend on the time span used, which can be reproduced 
with our model presented here (Zhang & Yi 2012). However a clear-cut and robust test 
of our model is to compare the prediction of the above equation. As discussed above, the 
first or the second term in the above equation may be ignored for old millisecond pulsars or 
young pulsars, respectively. The current on-going programs of intensive long term precise 
timing monitoring of pulsars should produce the required data for validating our model. In 
doing so, both the decay index and the oscillating properties of each individual pulsar can 
be determined, which in turn can provide important insights into the physics of NSs. 
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Fig. 3. — Correlation between — z> and v (panel (a)) and the exponential decay constant 
(panel (b)) for young pulsars with r c < 2 x 10 6 yr and v > 0. In panel (a), the dotted and 
dashed lines are the best-fit prediction of the exponential decay model, and the best fit with 
v oc — z> 7 (7 = 1.22). Panel (b) shows that a single distribution of t& is inconsistent with 
data. 
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Fig. 4. — Relations between a and rms of the fit (panel (a)), r] m (the median value of rj; 
(panel (b)), and t (the beginning time of power-law decay; panel (c)) for different initial 
values of B , for young pulsars with r c < 2 x 10 6 yr and i) > 0. The data requires a > 0.5. 
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Fig. 5. — Correlations between v and (— z>) 1+ ^/z/ 3/3 with a = 0.5 (panel (a)) and a = 1.0 
(panel (b)) for young pulsars with r c < 2 x 10 6 yr and v > 0. The dotted lines are the 
best-fit of v = 7](-u) 1+l3 /u^. 



-18 - 




Fig. 6. — Correlations between and i>z/ 3/3 for all pulsars with a = 0.5 (panels (a, b)) 
and a = 1.0 (panels (c, d)), in comparison with model predictions with different values 
of /. The solid lines correspond to / = 0. The left two panels show data with > 0, 
where the dotted and dashed lines are for = r}{— z>) 1+/3 /z/ 3/3 — 2z>/ (since v < 0) and 
v = v) 1+ P/v 313 + 20 f, respectively. The right two panels show data with < and 
model predictions of v — rj(— z>) 1+/3 /i/ 3 ^ + 20 f. The values of / are given in Figs. (jTJ) and 
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Fig. 7. — Correlations between / and — z> (panels (a, c)) and v (panels (b, d)), for a = 
(panel (a)) and a = 1.0, respectively. No significant correlation is found. 
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Fig. 8. — Distributions of f /a of all pulsars (black solid lines), pulsars with positive v 
dot-dashed lines) or negative v (green dashed lines), for a = 0.5 (panel (a)) and a = 
(panel (b)), respectively. The median values of f/a is about 3 x 10 -12 in both panels. 
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Fig. 9. — Comparisons between the observed v and that predicted by v = i](—i>) 1+l3 / h> 3/3 ±2i , f 
with f /a = 3 x 1CT 12 for a = 0.5 (left panel) and a = 1.0 (right panel). 
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Fig. 10. — Distributions of the braking index n b = vvjv 2 . The solid and dotted lines show 
n b > and n\, < 0, respectively. The vertical dashed line indicates nb = 3. The numbers of 
positive or negative values of n b are 193 and 173, respectively. The numbers of nb > 3 or 
n b < 3 are 192 and 174, respectively. 



